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Context. Massive Galactic clusters (> 1000 Mg) exhibit a clear correlation between cluster density, size and age and can be sorted 
in two categories, i.e. starburst and leaky clusters. The reason for the existance of two types of massive clusters is an open question. 
However, the answer is probably connected to a different formation histories of the two types. 

Aims. In this study we concentrate onleaky clusters only and investigate possible formation scenarios and gas expulsion phase. 
Methods. This is done by using existing observational data and numerical results of embedded cluster properties. 
Results. Assuming that a clear correlation between cluster density, size and age exists, it is shown that the density-radius development 
over time forembedded clusters can be approximated by p » 100 * r^'^Mepc"'. The consequences for the star formation process in 
leaky clusters are discussed and found to favour an inside-out star formation scenario with an initially low but later accelerated star 
formation rate. It is shown how the leaky clusters form in a unique sequential manner and that rapid gas expulsion is responsible for 
the 80-90% mass loss over the next 20 Myr. 
Conclusions. 
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C/^ 1. Introduction 
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Most stars form in dense clusters embedded within giant molec- 
C3 ular cloud clumps( |Lada & Lada||2003t|Pudritz||2002| i. Given the 
relatively high number of clusters younger than 10 Myr com- 
pared to the much lower number of older clusters in equivalent 
age spans, the majority of clusters must rapidly disolve soon af- 
ter their formation. As no clear correlation between cluster mass 
and age were found, until recently the general view has been that 
in the Milky Way, star formation results in a multitude of clus- 
ter with hardly any constraints on stellar content, size or density. 
However, at least for clusters more massive than 1000 M© this 
view had to be revised. [Pfalzner| ( |2009[ l showed that these mas- 
sive clusters follow two separate sequences. Although these are 
difficult to distinguish in the mass-age plane they become im- 
mediately apparent in density-radius space. In both cases the ex- 
pulsion of gas from the cluster via stellar winds, ionisation and 
supernovae explosions(Hills' '1980', 'Goodwin & Bastianl 2006 
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[Bastian et al. , 2008; Baumgardt and Kroupa 2007) leads to an 



expansion (see Fig. 1 here and more details in |Pfalzner ( 2009| l 
Fig. 2). 

The reason of this bi-modal development is currently not un- 
derstood. One might think that a difference in cluster mass might 
cause it. However, Fig. 1 shows that at an early age (<4Myr) both 
types of clusters have nearly the same mass but the cluster radii 
in the two groups (marked by red diamonds) differ by more than 
a factor 10. An additional difference is that the mass of starburst 
clusters remains more or less constant whereas that of leaky clus- 
ters reduces to about 10% - 20% of its initial value. A hint to the 
origin of the two groups of massive clusters is the obvious dif- 
ference of the location of the clusters - all starburst clusters are 



located in regions of high stellar density - near the galactic cen- 
ter or in the spiral arms - whereas the leaky clusters are situated 
in lower density areas. Clearly the answer to the question of why 
massive clusters exist in two categories must lie in their forma- 
tion process. In this paper we concentrate on the formation and 
development of leaky clusters which start out with masses of the 
order of 10"^ M© and loose 80%-90% of this mass over the next 
20 Myr after gas expulsion. 

Leaky clusters form as a bound entity which dissolves to 
a large degree after gas expulsion. The exphasis of this study 
lies on the formation process of leaky clusters. The precursors to 
leaky clusters must be found among the still embedded clusters. 
But the question is which embedded clusters develop into leaky 
clusters and how does this development proceed. 

One possibilty would be that leaky clusters form by merg- 
ing of sub-clusters like recently suggested by several authors for 
less massive cluster (see, for example, Kiipper et al. 201 1 and 
references therein). However, to form a leaky cluster this would 
require that several clusters like the Orion Nebula cluster (ONC) 
with > 1000 stars or hundreds of small-N clusters would have 
to form in close vicinity and merge within < IMyr. So far there 
exist no observational indications for such a formation process 
for leaky clusters. Here we stick to the more conservative picture 
of a leaky cluster forming as a single entity. 

In the following existing data of embedded cluster proper- 
ties are used to constrain the conditions of leaky cluster forma- 
tion and the gas expulsion phase. The idea is that a better under- 
standing of the formation process and later development of leaky 
clusters will eventually answer the open question why leaky and 
starburst clusters form and develop differently. 
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Fig. 1. Cluster mass as a function of cluster size (full symbols) 
for exposed clusters more massive than IO^Mq. The diamond 
symbols represent clusters younger than 4 Myr, the squares clus- 
ters in the age range 4Myr < < 10 Myr and the circles clusters 
with ages lOMyr < < 20 Myr. The thick drawn line shows that 
the cluster mass remains more or less constant for the starburst 
cluster sequence, whereas the leaky cluster mass rapidly declines 
with cluster a ge (thick dashed li ne) . The values were taken fr om 
Figer| poos) , [Wolff et al.| ( [2007l l and |Borissova et al.| ( p008] l. In 
addition the relation between the cluster density and cluster ra- 
dius are shown as in Pfalzner ( 2009 | l by empty symbols. 



2. Method 

This work uses available observational data of the properties of 
embedded cluster by Lada & Lada ( 2003| l to propose a possible 
formation history for leaky clusters. In principle one could de- 
duce the star formation history of massive clusters by looking at 
the age distribution within single such clusters. However, there 
are two problems with this approach: i) leaky clusters show con- 
siderable mass loss within the first 10 Myr of their developement 
and ii) it is not clear how long the process of star formation actu- 
ally takes in such clusters - it could be coeval or take up to ~10 
Myr 

In young embedded clusters like the ONC, wide spreads in 
colour-magnitude diagrams (CMD) have been observed, which 
have previously been interpreted as age spreads of up to 10 Myr 
in association members(Pall a and Stabler] [TOOO). However, other 
authors have suggested that the observed CMD spreads are due 
to a mixture of photometric errors, variability in pre-MS stars 



(Herbst et al. 1994 1, binarity (Preibisch and Zinnecker 1999 1, 



episodic accretion phases! Baraffe et al. 2009| or a spread in 



line-of-sight distances or extinctions. Currently age spreads of 
1-3 Myrs are regarded as realistic for the ONC. However, the 
age of Cyg OB2 was originally estimated as 1-3 Myr consider- 
ing only OB stars. A concentration of 5-7 Myr old A stars was 
found just to the south of Cyg OB2 (Drew et al. 2008) while 
an unclustered spread of evolved stars up to 10 Myrs old was 
identified over a wider area! Comeron et al. , 2008 ). Wright et al. 
( |2010( ) favour the view that the present-day OB stars are prod- 
ucts of the latest phase of star formation, but reason that at least 
one more epoch of star formation took place beforehand. In sum- 
mary, the observations of the stellar ages in single star clusters 
give a rather unclear picture of the duration and the dynamics of 
the star formation process. 



So rather than looking at the star formation history in one 
single cluster one could look at clusters of different age and con- 
struct a star formation history from this sequence. Mature leaky 
clusters consist of thousands of stars, but nevertheless they must 
have started out in the embedded phase containing just a few 
dozen stars, aquiring additional members via star formation. The 
straightforward approach would be to determine the ages of the 
embedded clusters and sort them temporally to obtain the for- 
mation history of leaky clusters. Unfortunately this is hindered 
by the fact that age determination in the embedded phase is even 
more complicated than in the early exposed phases due to the 
high uncertainty of the pre-main sequence (PMS) stellar evo- 
lution models. So although embedded clusters are usually at- 
tributed ages in the 0-2 Myr range, the error bar is generally of 
the same order. As a consequence ordering embedded clusters 
according to their age would be extremely error prone. 

For the exposed phase an important consequence of the se- 
quential nature of the leaky cluster development in the cluster 
radius-age space is that in this case the radius of a cluster can 
be used to determine its age (see Fig. 3 in 'Pfalzner ( 20091). For 
leaky clusters a fit to the observational data gives a correlation 
between the average cluster radius and its age td of the form: 



f,,/Myr ~ 0.15(rrf/pc) 
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(1) 



One way out of the dilemma of determining a formation his- 
tory for leaky clusters is the assumption that, similar to the ex- 
posed phase of the leaky cluster development, there exists a sim- 
ilar correlation between the cluster radius and its age in the em- 
bedded phase. In other words, we postulate a relation between 
cluster membership i.e. of the number of stars akeady formed, 
and age in the embedded phase for star clusters that eventually 
develop into leaky clusters. This is the simplest assumption we 
can make: it implies that star formation takes place in these clus- 
ters. 

Now implicitly the same underlying reasoning as for ex- 
posed clusters is used: Since each of the exposed leaky clus- 
ters has the same developmental history, it seems reasonable 
to assume that their precursors in the embedded phase have as 
well a common history. However, not every embedded cluster 
will develop into a massive leaky clusters. Only if there is suf- 
ficient matter in the surrounding cloud this developmental path 
is open. Clusters with too low mass content in their surrounding 
stop forming stars as soon as they run out of material - well be- 
fore they have reached masses in excess of 10^^ Mq. Most likely 
these starved clusters quickly disperse so that their surface den- 
sity more or less immediately drops below the detection line. 

Back to the formation of leaky clusters - in order to inves- 
tigate the consequences of the assumption of a common history 
of the precursors of the leaky clusters, we use the data of em- 
bedded cluster properties by Lada & Lada (2003). They selected 
these clusters within 2kpc from the sun in such a way that the 
group must have a stellar-mass volume density sufficiently large 
to render it stable against tidal disruption by the Galaxy (i.e., 
p >0. 1 Mopc"^) and by passing interstellar clouds (i.e., p >1 
Mopc"^). In addition, the group should contain a sufficient num- 
ber of members to ensure that its evaporation time is greater than 
the typical lifetime of Galactic open clusters. It means that pos- 
sibly existing clusters with densities < 5Mo(rc//pc)"^ are not 



included in this sample given in Lada & Lada ( 2003| l(see dashed 
line in Fig. 2). 

[Carpenter et al.| ( |2000| l state that for their sample the ob- 
servational density hmits imply that clusters of density < 31.5 
M0pc"^(rc;/pc)"' are too tenuous to be detected against the 
background of the field. If this limit also applies to the |Lada &| 
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Fig. 2. Cluster density as a function of cluster size for the embed- 
ded clusters listed in Lada & Lada| f2003 ) (coloured symbols) 
and the leaky clusters from Pfalzner 2009 (black symbols). The 
open coloured symbols for the embedded clusters indicate the 
values for IRAS 06155+2319, GEMl, NGC 2281 and GEM4. 



The selection criteria by Lada & Lada (2003) imply that clus- 
ters with densities < 5Mo(rci / pcy~^ are not included. In addition 
clusters with densities < 31. 5MQpc^^(rci / pcy^ can not be de- 
tected against the field of background stars. 
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Fig. 3. Mass as a function of radius for embedded clusters. The 
data points are taken from [Lada & Lada (2003 1 and the best fit 
for a/?* gives 0=359+15 and i>= 1.7 1+0.07. If the cluster mass 
is corrected to a lower mean stellar mass, this would change to 
flT?* gives a=451±15 and i?=1.71±0.07. However, given the large 
errors a/?^-relation could as well be considered. The symbols are 
the same as depicted in Fig. 2. 



|Lada| ( p003] l data sample (see dashed line in Fig. 3.), the absence 
of lower density embedded clusters might be an artifact of obser- 
vational limitations. On the other hand, this limit might also be 
too stringent - the exposed leaky clusters with sizes >10pc can 
be detected despite their density being relatively low compared 
to the background density. 



ably and the cluster radius nearly grows by a factor of 10. During 
this phase the cluster density stays nearly constant or decreases 
only slightly. Whether these are actually two distinct phases or 
a continuous transition cannot be decided based on current data. 
However, considerable acceleration of star formation occurs at a 
cluster size of about 0.7- Ipc when more than 100 members are 
formed. In Fig. 2 the dependency of the cluster density p on the 
cluster radius r^; 



3. Results 

3. 1 . Density development in fhe embedded pliase 

Even the most massive leaky cluster must have started out at 
some point with a small number of stars. Fig. 2 shows the lower 
part of Fig. 1 - the leaky cluster sequence - but augmented by 
the data sample of embedded clusters from |Lada & Lada ( 2003| l 
with known sizes and masses. Because this data set is limited 
to embedded clusters within 2 kpc from the sun, its members 
are sufficiently far away from the Galactic center and the spiral 
arms, so it is most likely free from precursors to starburst clus- 
ters. Fig. 2 implies that the precursors of leaky clusters form an 
extension of the leaky cluster sequence towards smaller radii (< 
3pc) at density of > 100 Mopc However, the embedded clus- 
ters do not form a direct continuation of the slope of the main 
evolutionary sequence to lower radii, but the density develop- 
ment is flatter. This is because expansion due to gas expulsion 
has not started yet and star formation is still ongoing. 

In Fig. 2 and Fig.4 clusters with <50, 50-100, 100-300 
and >300 members have been assigned to four different groups 
via distinct colours. The average radius of these cluster groups 
increases with the number of members. Following the time- 
sequence hypothesis - that age and radius are also correlated 
in the embedded phase - we may interpret Fig. 2 as follows : 
A cluster starts with an initial phase where a limited number of 
stars form (< 100), the cluster radius increases by a factor 2-3 and 
the cluster becomes more diffuse. The relatively high number of 
embedded clusters is indicative of the slow star formation rate at 
this stage. Afterwards the star formation rate increases consider- 



lOOr: 



(2) 



provides only a rough clue to the development because the 
in the original data set not specified but presumably large errors 
in cluster radius do not justify a proper fit at this point. The error 
in mass has been estimated as approximately a factor 2 by Lada 
& Lada (2003), this error was taken into account in Fig. 2 for 
determining the stellar densities. 

Fig. 2 shows that the average density in low N-clusters - 
which in the here presented picture corresponds to the youngest 
clusters - is often more than 10 times higher than that of the 
ONC. This means that in these early phases nearly all stars 
are exposed to interactions with neighbouring young forming 
stars which most likely influences star and planet formation pro- 
cesses. Properties that might be altered by this high stellar den- 
sity are, for example, the binary frequency, the density profile in 
the circumstellar discs, etc. Although clusters like the ONC show 
a similar density in the central cluster regions to that of low-A^ 
clusters, most stars are located in areas where the environment 
is less important. 

The relative importance of these environmental effects on 
star and planet formation is still an open question, however, the 
here provided knowledge of the density development during the 
star formation process allows to determine the encounter proba- 
bility in such clusters better in future investigations. In the future 
this can show in how far the high stellar density influences star 
and planet formation processes in leaky cluster environments. 

As mentioned before, not all embedded clusters included in 
Fig. 2 will necessarily develop into leaky clusters: Although the 
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|Lada & Lada| ( |2003[ ) sample probably excludes precursors to 
starburst clusters, it does not exclude precursors of lower mass 
clusters. Only if there is sufficient material in the surrounding 
cloud will the cluster continue to evolve towards the leaky cluster 
sequence. It is, for example, unlikely that IC348 at the leftward 
end of the embedded track in Fig. 2 will develop into a leaky 
cluster because its gas content is too low. For embedded clusters 
that develop into lower-mass exposed clusters the relation of Eq. 
2 probably serves as an upper density limit. 

However, clusters that contain sufficient gas will eventually 
join up with the leaky cluster sequence when the gas expulsion 
has caused an expansion to ~ 2-4pc. It is most likely that only 
embedded clusters as dense as those shown in Fig. 2 can be the 
percursors of the < 4 Myr old leaky clusters. Lower density clus- 
ters might exist, but those that are precursors to leaky clusters 
will be at the high density limit (for a given radius) of this data 
set. 

In Fig. 2 the dependence of the cluster density on the cluster 
radius could as well be interpreted as an observational selection 
artifact, in the sense that many clusters show multiple peaks in 
their density. So if the low-mass clusters would all represent one 
such core and the higher mass clusters several cores with low 
density regions in between, one could obtain a similar diagram. 
However, Fig. 3 shows clearly that it is really a mass-radius rela- 
tion. 

Plotting the cluster mass of the Lada & Lada ( 2003[ l data as 
a function of the cluster radius (see Fig. 3) a relation of the form 



(3) 



with fl=359+15 and Z7=1.71±0.07 is obtained. It can be seen that 
the extension of Eq. 3 to larger cluster radii or equivalent in later 
times in the cluster development, meets with the leaky cluster 
track at the point where these clusters become exposed. So Fig. 
3 represents the formation, gas expulsion and dissolution phase 
of the leaky cluster sequence. 

In the following a rough estimate of the errors in cluster radii 
is attempted. There are several sources that could lead to such 
an error. Due to the embedded nature not all stars are detected. 
If mass segregation is primordial, the higher likelihood of low- 
mass stars to be missed could lead to a too small detected radius. 
In addition, many of the embedded cluster are not spherically 
symmetric but are more extended in one direction. In Fig. 3 it 
was assumed that the error in cluster radius is approximately 
half its actual value. 

[Adams et al.| ( |2006l l noted a related correlation between clus- 
ter radius and number of stars A^, namely, rd tx:= AT'^^. However, 
they interpret this correlation as a static property of embedded 
clusters and not as an outcome of temporal development of the 
cluster formation process. Nevertheless, with M oc N it would 
follow that Mc oc which considering the errors could as well 
fit the data (see Fig. 3). 

If the trend expressed by Eq. 3 is not an evolutionary se- 
quence at all, then one has to explain why only certain mass- 
radius combinations seem to exist for young embedded clusters. 
Selection effects could play a role, so it would be hard to detect 
clusters with Mc < 100 and r^i < 2pc. However, it would 
be unlikely to miss clusters with, for example, Mc > 500 Mq 
and rd ~ 0.5pc, which are clearly not present in the sample. One 
way out, would be a mass-size relation for cores (see, for exam- 
ple. Swift & Williams 2008). However, then the percursors of the 
all clusters would have to form from thinly spread clusters with 
very low initial surface densities and be all below the detection 
limit. These clusters would then form stars within a more or less 



fixed radius finally reaching such a density that they can be de- 
tected. This explanation works very well for low-mass clusters. 
However, in such a picture the precursors of the leaky clusters 
should be embedded clusters with cloud masses Mdoud > 1000 
Mq containg thousands of stars within rd ~ 3-6pc. Such clusters 
should be detectable in the solar neighborhood. As this is not the 
case, the only explanation would be that they are either rare or 
very short-lived. 

Therefore we prefer here the model of interpreting the mass- 
radius relation as a time-sequence. In this model the existence of 
a mass-radius relation naturally leads to such a form as described 
by Eqs. 2 and 3 because these clusters still generate stars, so log- 
ically those with smaller populations must simply be younger. 

Eq. 3 probably directly describes the star formation process 
in leaky clusters since loss processes such as escapers play only 
a minor role. ' Weidner et al.| ( |201 1) showed that in the embed- 
ded phase only about 15-20% of the cluster mass is lost through 
escapers with little consequence to the cluster radius. If a direct 
relation between cluster radius and cluster age exists, Eq. 3 im- 
plies strongly accelerated star formation. 

3.2. Star formation history 

What does the embedded cluster sequence tell us about the star 
formation in a cluster environment? Two star formation scenar- 
ios would be possible: i) either stars form only in the central 
high-density area, or ii) the star formation process itself pro- 
gresses from the inside outwards. We assume that the cluster 
size is a function of time of the form rd where c is a con- 
stant with 0< c < 1 . Limiting our investigation to c < 1 means 
an exclusion of acceleration cluster expansion. Eq. 3 translates 
to Mc oc f'^'^. In the case i) the velocities would have to be su- 
pervirial so that as time progresses the cluster size expands. This 
could possibly be caused by an early gas expansion. However, 
this would require all clusters to have the same velocity disper- 
sion. 

Case ii) is easier to reconcile with above findings. In a ho- 
mogenoeus medium, one would expect Mc oc if star formation 
just proceeded linearly from the centre outward. Since the ob- 
servations indicate a Mc oc f '^'-dependence, this would either 
require the star formation front to decelerate with time rd oc f"-^'" 
or, the medium to become less dense. The latter fits what we ex- 
pect from theory as well as observations - the gas density in the 
cluster forming clumps decreases from the centre outward. 

The gas densi ty pg = M{r) /V in young forming clus- 
ters! Larson 1981 Chandrasekhar 1939| l is either proportional 
to r"^ or r"'^. Besides, power law density profiles with slopes 
in the range 1.5-2.0 are put forward for molecular clumps by 
various studies (see |Parmentier| pO ll\ and references therein). 
This means that the matter available for star formation has the 
same dependency. In the first case one would end up with the 
gas mass increasing as Mc oc and the second case would result 
in a dependency Mc oc r'. Assuming a constant star formation 
efficiency throughout the cluster and using the size-age relation, 
this corresponds to the cluster mass increasing as Mc oc t^'^ and 
Mc oc respectively. Generally, one could conclude that for a 
linear cluster expansion (c=l), the gas density distribution in the 
clump from which these clusters form would have to be of the 
form Mc oc r'-^. 

An alternative picture would be that it is not the density 
distribution that is responsible for the shape of the observed 
density-age relation, but the way that star formation proceeds. 
Here it is assumed that star formation continues in the central 
area while the star formation front moves outwards in a medium 
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cluster age [Myr] 

Fig. 4. Cluster size development as a function of cluster age 
for exposed clusters. The observational values were taken from 
Wolff et al. (2007). The simulation data points are from Goodwin 
& Bastian (2006). a linear fit in the log-log space is shown as 
drawn line, scaled versions of the results of Goodwin & Bastian 
(2006) for instantaneous gas expulsion with 10% (short dashed 
line), 30%(long dashed line) and 50%(dotted line) star formation 
(SF) efficiency are shown. 



with a r"^ gas density dependence. This is a likely scenario, as 
star formation would probably not immediately stop in the cen- 
tre but continue until enough massive stars have formed to blow 
a cavity into the gas in the central area. For simplicity in Fig. 3, 
a model assuming linearly propagating star formation has been 
calculated (dashed line). In this picture in each shell of thickness 
Ar initially the same amount of star formation AM happens, but 
the further a shell is from the center the later it starts to produce 
stars. After its onset the star formation rate slows down by the 
factor ar in each of the individual shells as the gas content in 
these regions decreases. Assuming that the star formation pro- 
ceeds over N shells, the total star formation Ms f is given by 



Msf{N) 



N N 
i=i 1=1 



i+ l)a'f"'+'AM. 



(4) 



This model fits well with the derived density radius relation 
shown in Fig. 3 if we choose arit„+i) = 0.7 * ar{t„). 

These simple estimates should be regarded as a first attempt 
to show how new insights on embedded cluster dynamics might 
provide clues to the star formation process itself This should 
motivate future observations as well as more comprehensive nu- 
merical modelling. 

3.3. Gas expulsion 

As Fig. 2 demonstrates there are two distinct phases in the den- 
sity development of leaky clusters - the star formation phase 
which leads to a density/radius relation as expressed by Eq. 2 
and the cluster expansion phase which is described by Eq. 1. We 
stress that while the sequence of leaky clusters stems from dy- 
namical evolution (i.e. post-gas expulsion expansion), the em- 
bedded cluster track is the imprint of the formation process. The 
"switch" from one to the other process is the gas expulsion as 
illustrated by Fig. 3. 



Since the cluster expansion and mass loss are the result of 
the star formation process and the subsequent gas expulsion, it 
should be possible to draw conclusions about the early history of 
these clusters from the properties of these two processes. Using 
results from existing numerical simulations we can place con- 
straints on how the gas expulsion is likely to proceed in leaky 
clusters. 



Goodwin & Bastian ( 2006 1 happened to simulate a situation 



similar to that in leaky clusters. Their cluster simulations con- 
tained 15 000 stars distributed according to a Plummer modell 
with a core radius of 1 pc. They investigated the process of in- 
stantaneous gas expulsion by assuming virial equilibrium before 
gas expulsion and varied the star formation (SF) efficiency. Here 
SF efficiency is defined as M stars /(M stars + Mgas) with M stars be- 
ing the mass of the stars in a cluster and Mgas the gas mass. Note, 
that leaky clusters have a somewhat larger number of stars, typ- 
ically ~ 20000-40000. 

Fig. 4 shows the radial development of leaky clusters as a 
function of cluster age for the data used in Pfalzner ( 2009| l. Here 
the drawn line represents Eq. 1 which is basically a linear fit in 
the log-log space. Note the relatively large errors of 1-2 Myr in 
cluster age. In addition Fig. 4 shows the simulation results of 



Goodwin & Bastian ( 2006/ scaled to the initial cluster param- 
eters relevant here. Such a scaling is possible since the model 
excludes scattering processes between the cluster stars. It can be 
seen that a low SF efficiency of 10% leads to somewhat smaller 
cluster radii than observed whereas a high SF efficiency of 50% 
would lead only to an increase of cluster radii by a factor of ~ 2. 
The best fit is obtained by an instantaneous gas expulsion with a 
30% SF efficiency. In this case the simulated and observed leaky 
cluster expansion are consistent. Note that the data of the 30% 
SF efficiency simulation actually seem a better fit to the observa- 
tional data than Eq. 1 . However, the age determination has such 
a large error that it is not possible to decide which one is more 
realistic. Only better observational data and simulations tailored 
to the actual parameters of leaky clusters will give an answer to 
this question. 

It has been suggested that the distinction between starburst 
and leaky clusters is basically identical to that of systems re- 
maining bound and becoming unbound after gas expulsion (| 
Portegies Zwart et al. 2010 Gieles & Portegies Zwart 201 l[ l 
and that the leaky clusters are therefore (unbound) associations. 
However, taking the above derived SF efficiency of 30%, the 
results oflGoodwin & Bastian| (|20061) indicate that despite the 
considerable mass loss of these clusters, about 10%-15% of the 
stars can be expected to remain bound and develop into a cluster 
which lasts > lOOMyr. Although much reduced in mass, these 
bound entities still consist of ~ 1000 stars albeit with a low vol- 
ume density of the order of 1 Mapc"^. It is a matter of cluster 
definition, if such a group of star is still called a cluster If one re- 
gards boundedness of a sufficient number of stars as prime prop- 
erty of a cluster then it is. By contrast, if one requires a certain 
volume or surface density (for a review of the different defini- 
tions see jBressert et al.| ( |2010| l) then it is a cluster at the start of 
the expansion but can no longer regarded as such at later stages. 
In either case more detailed investiagtions should be carried out 
to specifiy the properties of these remaining entities more pre- 
cisely. 

4. Conclusion 

In this paper we investigated the formation and early develop- 
ment of leaky clusters. This first approach is far from definitive, 
however, it does provide a consistent picture. We showed that the 
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difference between leaky and star burst clusters is not mass re- 
lated but must lie in their early development (see Fig. 1). Taking 
our cue from the size/age relation demonstrated for young leaky 
exposed clusters, we postulated an analogous relation for the em- 
bedded phase. Based on this concept, 

- we examine the star formation history of leaky clusters by 
determining the cluster radius development as opposed to 
the traditional examination of the stellar age distribution in 
single clusters. Here the global properties of a multitude of 
known embedded clusters is used to make deductions about 
the star formation history inside individual star clusters. This 
is only possible because recent findings strongly suggest that 
there exists a common developmental history for leaky clus- 
ters. 

- the radial cluster development implies that star formation 
initially proceeds at a low rate and accelerates considerably 
later on. Spatially this is consistent with star formation in a 
sphere of increasing volume. 

- the increase of cluster mass with cluster size follows a 
dependence of the form - ar'^^j with fl=359±15 and 
^7=1.71+0.07. and the cluster density develops as p » 
100r;/l 

- it has been shown that the observational results of cluster ex- 
pansion are consistent with a 30% star formation efficiency 
followed by an instantaneous gas expulsion. 

From the existing data, it is not absolutely clear, whether the 
star formation process is the result of two distinct generations of 
stars, or a continuous but accelerated star formation. From the 
mass-radius diagram (Fig. 3) the latter seems more likely. 

With these results a picture of the mass and density develop- 
ment of leaky clusters can be constructed over the first 20 Myr of 
their existance: Leaky clusters start out with several tens of stars 
being confined within ~ 0.2 - 0.375c, grow in mass and cluster 
membership as Mc ~ 360r'-^ 
mass of several 10"* M, 
sion sets in and the clusters expand rapidly while loosing most 
of their mass over the follwing 20 Myr Throughout the cluster 
formation phase the density (at least in the central region) is so 
high that the cluster environment might significantly influence 
the planet formation process. 

In this paper the focus was on the early evolution of leaky 
clusters, it would be highly desirable to obtain the equivalent 
information for starburst clusters. However, at present it is im- 
possible to determine in which way star formation in starburst 
clusters differs from that in leaky clusters. Fig. 1 implies that 
the precursors of starburst clusters would form an extension of 
the starburst cluster sequence towards smaller radii (< O.lpc) at 
density of > 10^"^ Mopc Otherwise it would require them to 
contract and increase in density while simultanuously forming a 
high number of massive stars. No massive embedded cluster has 
been definitely classified as a starburst precursor so far possibly 
with the exc eption of W43 ( Motte et a l. , 2003 ; Nguy en Lu'O'Ng 
201 1\ This is not surprising as starburst clusters are situ- 



^1 until they reach their maximum 
©with a size of ~ 3-5 pc. Then gas expul- 
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ated either close to the Galactic centre (8kpc) and/or along spiral 
arms. At these locations the higher background density makes 
cluster detection significantly more difficult, so that even most 
of the exposed massive starburst clusters were only detected re- 
cently. 
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